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Abstract 

A  simple  method  of  computing  long-wave  radiative  cooling  in  the 


M 


troposphere  associated  with  water  vapor  is  described.  Radiation  from 
ozone  and  carbon  dioxide  is  not  considered.  However,  influences  of 
arbitrary  vertical  distributions  of  cloud  and  moisture  are  included. 

Average  annual  cooling  rates  along  a  meridional  cross-section  are 
calculated  for  a  cloudless  atmosphere.  The  results  agree  fairly  well 
with  the  total  radiative  cooling  (long-  and  short-wave)  as  given  by 
Manabe  and  Moller  (1961)  except  in  the  lower  troposphere  at  low 
latitudes.  Here  short-wave  absorption  by  water  vapor  is  appreciable. 

Long -wave  radiative  cooling  is  also  computed  in  a  case  of  a 
developing  cyclone  for  comparison  with  release  of  latent  heat.  The 
largest  cooling  occurs  at  cloud  top  and  can  be  a  significant  fraction  of 
the  amount  of  energy  released  as  latent  heat  in  the  upper  troposphere. 
Destabilization  of  the  cloud  mass  and  subsequent  increase  in  precipi¬ 
tation  may  be  important  in  cyclone  development. 


1 .  Introduction 

The  heat  budget  of  the  atmosphere  may  be  partitioned  as  follows: 

(1)  absorption  of  short-wave  radiation,  (2)  sensible  heat  transfer  from 
the  earth's  surface  by  conduction  and  convection,  (3)  release  of  latent 
heat,  and  (4)  absorption  and  emission  of  long-wave  radiation.  Of  these, 
long-wave  radiation  is  the  most  important  term  in  the  mean  tropospheric 
heat  budget  (see,  e.g.,  Davis  1960  or  London  1957).  Release  of 
latent  heat  is ,  on  the  average ,  the  second  most  important  term .  How¬ 
ever,  it  is  often  locally  predominant  in  areas  of  precipitation  and  its 
inclusion  in  a  prediction  model  appears  desirable  in  forecasting  cy¬ 
clogenesis  (Danard  1964,  1966a,  1966b). 

This  paper  describes  a  simple  method  of  including  long-wave 
radiation  from  water  vapor  in  a  tropospheric  numerical  model  which 
predicts  the  field  of  moisture .  Radiation  from  ozone  and  carbon 
dioxide  is  not  considered.  However,  effects  of  arbitrary  vertical 
distributions  of  cloud  and  moisture  are  included. 

Using  a  numerical  model  incorporating  release  of  latent  heat  but 
not  long-wave  radiation,  the  author  (L'anard  1966a,  1966b)  noted  a 
tendency  for  predicted  heights  to  be  too  high  in  the  upper  troposphere. 
This  may  have  been  due  to  neglect  of  long-wave  cooling  from  cloud 
tops .  Inclusion  of  long-wave  radiation  should  extend  the  viability 
of  such  models . 


2.  Theory  and  basic  procedure 

Since  flux  divergence  associated  with  carbon  dioxide  is  usually 
small  in  the  troposphere  (see,  e.g.,  Manabe  and  Moller  (1961),  or 
Manabe  and  Strickler  (1964)),  only  water  vapor  will  be  considered. 
Computations  will  be  performed  for  the  atmosphere  below  some  level 
Pu  (near  the  tropopause) .  This  is  the  highest  level  for  which  meteor¬ 
ological  variables  are  assumed  known. 

Results  of  Kuhn  (1963c)  suggest  the  following  expression  for  the 
corrected  precipitable  wa':er  or  path  length  between  p=0  and  p=p: 


Here  p  and  T  oenote  standard  pressure  and  temperature,  and  q  is 
o  o 

the  specific  humidity,  w  may  be  regarded  as  a  vertical  coordinate 
increasing  downward. 

In  computing  w^,  the  value  of  w  at  p=p^,  use  is  made  of  Manabe 

and  Moller1  s  (1961)  observation  that  the  frost  point  in  the  lower 

stratosphere  is  about  190  K  irrespective  of  season  and  latitude.  This 

value  is  attained  at  a  height  of  about  15  km  or  120  mb.  Assume  that 

the  layer  above  p^  is  isothermal  with  temperature  T  .  Setting 

T  =220  K,  (1)  then  gives  a  value  for  w  from  p=0  to  p=120  mb  of 
U  -7  -2 

4.6x10  gm  cm  .  Assume  that  q  varies  linearly  with  pressure  from 
p=120  mb  to  P=PU-  With  the  aid  of  (1),  this  permits  us  to  compute  the 
contribution  to  w  from  the  layer  between  these  levels  .  Adding  this  to 
the  contribution  from  the  layer  above  yields  (in  cgs  units) 

w  =  4.7  x  10"5  +  9.9  q  (2) 

u  u 
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where  q  is  the  value  of  q  at  p=p  .  Note  that  w  normally  decreases 
u  u  u 

from  equator  to  pole . 

The  net  flux  (positive  downward)  at  a  level  p*  between  cloud 
layers  (see  Fig.  1)  is 


w 


w*  c 

F*  =  F.  -F.  +;  {F  -F.  k'(w*-w)dw+J  {F.  -F.  }t'(w-w*)dw 


ba  be  J  bw  ba  J  .  be  bw 

w  w* 

a 


(3) 


where  F  ,  and  F^r  are  the  black  body  fluxes  corresponding  to  the 
temperature  at  w,  w  and  w  ,  and  <'(w*-w)  is  the  rate  of  change  of 

d  C 

emissivity  c  with  path  length  evaluated  at  a  path  length  (w*-w) 

(corresponding  to  the  slab  between  w  and  w*).  Throughout  this  paper, 

arguments  of  functions  are  enclosed  in  parentheses  ( ).  A  similar 

equation  is  given  by  Kuhn  (1963a).  If  the  level  p*  is  in  cloud,  F*  is 

assumed  zero.  If  no  cloud  exists  above  p*,  w  =0  and  F,  =0.  Eq.  (3) 

a  ba 

is  derived  in  a  manner  similar  to  that  described  by,  e.g. ,  Haltiner 

and  Martin  (1957,  pp.  85-86).  An  outline  of  the  derivation  is  given 

in  the  remainder  of  this  paragraph.  The  flix  in  the  spectral  interval 

X  to  X  +  d  X  at  w*  from  the  black  body  surface  at  w  is 

a 


dF**  ■2rteXH3<kXtw*-w.J) 

3  A 


(4) 


00 


where  H  (x)  =  J  e  Xyy  ndy,  F,  .  Is  the  black  body  flux  at  w  in 
n  Dei  A  3 

o 

the  spectral  interval  X  to  X  +  d  X  ,  k  is  the  absorption  coefficient 
and  X  is  the  wave-length.  The  flux  at  w*  from  the  stratum  between 
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w  and  w*  is 

d 


w* 

dFwX=2kX^w  PbwXH2(kXtw*-wi)dW  (5) 


where  is  the  black  body  flux  at  w  in  the  spectral  interval  X  to 

X  +dX  .  Now 


w* 

H3(kx{w*-wal)  =  £-kx  J  H2(k^  {w*-w})dw  (6) 

Wa 

Adding  (4)  and  (5),  making  use  of  (6),  and  integrating  over  all  wave¬ 
lengths  yields 

w*  « 

F‘*  =Fte-2j  J  kxFbaXH2(kx[w.-w})dXdw 

w  o 
a 


w*  <® 

+  2  J  J  kxFbwXH2(kX^W*”W^dXdW  (7) 

Wa  ° 


Now  for  an  isothermal  slab, 

00 

c'(w*-w)  =  J  kxFbwXH2(kX^W*'w})dX  (8) 

Eq.  (8)  is  derived  by,  e.g.,  Brooks  (1950).  It  may  also  be  obtained 
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from  (5).  For  a  fixed  value  of  (w*-w),  <  and  therefore  €'  as  well  will 
be  assumed  independent  of  temperature.  This  is  because  the  temperature 
dependence  has  already  been  included  in  (1),  which  is  used  to  compute 
w  and  w*.  Hence  (8)  will  be  assumed  to  apply  also  when  the  temp¬ 
erature  varies  between  w  and  w*.  Thus  the  last  term  on  the  right  side 
w* 

of  (7)  is  p  p  c.(w*_w)dw  •  Making  the  approximation 
rbw*  v 

W 

a 


F,  ./F.  =  F,  ./F,  (i.e. ,  for  a  fixed  wavelength,  F../F,  is 

ba  X  ba  bw  X  bw  b  X  b 

assumed  to  vary  little  over  the  range  of  tropospheric  temperatures),  the 

w* 

second  term  on  the  right  side  of  (7)  becomes  p  j  fi(w*  w)dw* 

w 

a 

Hence  (7)  gives  the  first  and  third  terms  on  the  right  side  of  (3) . 

Similar  calculations  show  that  the  contributions  from  the  black  body 
surface  at  w  and  the  layer  between  w*  and  w  are  represented  by  the 

C  o 

second  and  final  terms  on  the  right  of  (3) . 

The  assumption  that  a  cloud  surface  behaves  as  a  black  body  is 
probably  justified  for  water  droplet  clouds  at  low  elevations.  For 
example,  Kuhn  (1963a)  measured  the  emissivity  of  an  overcast  strato- 
cumulus  cloud  deck  with  base  at  900  mb  and  top  at  855  mb.  He  ob¬ 
tained  the  value  of  0.98.  As  Kuhn's  results  indicate,  however,  the 
attumption  is  less  realistic  for  clouds  at  higher  elevations.  For  a 
9/10  sky  cover  of  Jirrostratus  with  base  at  235  mb  and  top  at  160  mb, 
the  computed  emissivity  was  0.59. 

The  heating  rate  associated  with  radiation  is  given  by 


=  _  2_ 
v  at  J  c 


9  F* 
9p 


(9) 
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3 .  Adaptation  to  a  numerical  model 

For  purposes  of  illustration,  the  model  described  in  Section  .1  will 
be  adapted  to  a  primitive  equations  numerical  model  developed  by 
Krishnamurti.  However,  it  could  be  modified  to  any  model  in  which  e. 
moisture  parameter  is  predicted. 

The  level  p^  is  taken  as  200  mb  and  Eq.  (2)  is  used  to  compute 
w^.  Eq.  (3)  is  evaluated  at  200,  400,  600,  800  and  1000  mb  (Fig.  2). 
The  integrals  are  computed  by  summation  over  200  mb  thick  layers. 

For  example , 

w* 

I  tFbw"Fbai€,(w*'w)dw  =  E  tFbw-FbaHc(w*-wt)-c(w*-wb)}  (10) 
w 

a 

Here  w^  and  w,  are  the  values  of  w  at  the  top  and  bottom  of  the  200  mb 
t  b 

thick  contributing  layer,  and  Fbw  is  the  black  body  flux  at  the  mid¬ 
point.  The  summation  extends  over  all  layers  between  w  and  w*. 

d 

Similarly, 

w 

J  C(Fbc-Fbw}€'(w-W*)dw  =  S  CFb  -Fbw)Cc(wb-w*)-c(w  -w*)}  (11) 

w* 

In  (1 1) ,  the  summation  extends  over  layers  between  w*  and  w  .  Kuhn 

c 

(1963b)  employed  a  similar  method  of  computation. 

Cloud  is  assumed  to  exist  at  a  given  level  if 

q/q  >r(p)  (12) 

where  q  is  the  saturation  specific  humidity  and  r(p)  is  a  critical  value 
s 

for  existence  of  cloud.  Table  1  gives  values  of  r(p)  adapted  from 
Smagorinsky's  (1960)  relations  between  reladve  humidity  and  cloud 
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Table  1 .  Values  of  r(p)  (see  Eq .  (12))  for  existence  of  cloud  (adapted 
from  Smagorinsky  1960). 


p(mb)  900  700  500  300 

r(p)  0.84  0.75  0.70  0.69 


[ 

I 

amount.  The  values  in  the  table  correspond  to  8/10  cloud  at  the  given 
level.  Relation  (12)  is  tested  at  300,  500,  700  and  900  mb,  with  the 
(  cloud  (if  it  exists)  corresponding  to  each  level  assumed  to  be  200  mb 

thick.  The  model  atmosphere  can  thus  have  sixteen  configurations  of 
cloud  ranging  from  completely  cloud- free  to  cloudy  at  all  four  levels 
(Fig.  3) 

Kuhn's  (1963c)  values  of  c  are  used. 

In  the  computer  programming  of  Eq.  (3),  F*  is  first  set  equal  to 
zero  at  the  five  levels  indicated  in  Fig.  2.  Then  it  is  evaluated  from 
200  mb  down  to  the  top  of  the  highest  cloud  (or  the  earth's  surface). 
This  means  that  at  this  stage  flux  computations  are  complete  for  con¬ 
figurations  1-5  of  Fig.  3.  Next,  F*  is  computed  from  the  base  of  the 
highest  cloud  down  to  the  top  of  the  next  highest  cloud  (or  the  earth's 
surface).  This  completes  calculations  for  configurations  6-15. 

Finally,  F*  is  computed  at  the  800  and  1000  mb  levels  of  configuration 
16  to  finish  the  flux  computations. 

Eq.  (9)  is  evaluated  at  300.  500,  700  and  900  mb  by  expressing 
5F*/dp  as  a  finite  difference  over  layers  200  mb  thick  (see  Fig.  2). 
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4.  Application  to  a  mean  meridional  cross-section 

Manabe  and  Moller  (1961)  compute  the  annual  average  cooling 
rates  associated  with  radiation  as  a  function  of  height  and  latitude 
for  the  northern  hemisphere.  They  include  long  and  short-wave  ab¬ 
sorption  by  ozone,  carbon  dioxide  and  water  vapor  and  neglect  the 
effects  of  clouds. 

Since  only  long-wave  absorption  by  water  vapor  is  included  in 
the  method  presented  in  this  paper,  it  is  of  interest  to  compare  the 
tropospheric  cooling  rates  thus  computed  with  those  of  Manabe  and 
Moller. 

Using  the  mean  values  of  temperature  and  pressure  compiled  by 
London  (1957),  Eq.  (3)  was  evaluated  at  1  km  intervals  from  the  sur¬ 
face  to  12  km.  London's  humidity  values  were  accepted  for  T  >  -30  C. 
However,  at  lower  temperatures,  a  constant  frost-point  lapse  rate  of 
6.25  C  km  1  was  assumed.  This  was  also  done  by  Manabe  and  Moller. 

Eq.  (2)  was  used  to  compute  w  and  p  was  taken  as  200  mb.  For 

u  u 

comparison  with  Manabe  and  Moller's  results,  the  atmosphere  was 
assumed  cloud-free.  Cooling  rates  were  computed  at  0.5,  1.5,...., 
11.5  km.  Average  annual  values,  which  have  been  smoothed  slightly, 
are  presented  in  Fig.  4.  Fig.  4  may  be  compared  with  Fig.  21  of 
Manabe  and  Moller  (1961).  The  two  figures  are  fairly  similar.  The 
average  cooling  rate  in  Fig.  4  is  1.1  C  day  *  compared  to  0.9  C  day  * 
for  Manabe  and  Moller's  results  for  the  same  region.  Differences 
between  the  results  of  this  paper  and  those  of  Manabe  and  Moller  may 
be  seen  in  Fig.  5.  It  will  be  seen  that  considering  just  long-wave 
radiation  from  water  vapor  tends  to  overestimate  the  total  radiative 
cooling  (long-  and  short-wave)  at  low  latitudes  and  low  levels.  The 
discrepancy  is  probably  due  to  short-wave  absorption  by  water  vapor. 
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Results  of  Davis  (1963)  indicate  that  in  such  regions  this  absorption 
can  be  half  as  large  as  the  long- wave  cooling. 

Goody  (1964,  p.  272)  states  that  cooling  rates  computed  from 
different  techniques  may  vary  by  as  much  as  50  per  cent.  In  view  of 
this  remark  and  the  assumptions  made  in  the  present  paper,  the 
agreement  with  Manabe  and  Moller’s  results  must  be  considered  as 
reasonably  good.  Note  that  the  above  comparison  is  not  simply  a 
test  of  the  method  used  to  compute  long-wave  radiation  from  water 
vapor.  It  is  primarily  an  assessment  of  the  validity  of  omitting  short¬ 
wave  radiation  as  well  as  long -wave  radiation  from  carbon  dioxide 
and  ozone  in  the  troposphere. 

Fig.  4  may  also  be  compared  with  the  infrared  cooling  rates  of 
London  (1957).  However,  London  accounted  for  average  cloudiness 
whereas  in  the  present  study  the  atmosphere  was  assumed  cloud-free. 
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5  .  Comparison  with  release  of  latent  heat 

The  author  has  previously  calculated  the  spatial  distribution  of 

release  of  latent  heat  at  400,  600  and  800  mb  in  a  developing  cyclone 

over  the  United  States  (Danard  1964).  To  permit  a  comparison  between 

these  results  and  long-wave  radiative  effects,  the  method  described  in 

Section  3  was  modified  somewhat.  The  level  p  was  taken  as  300  mb 

u 

and  T^  was  assigned  the  value  230  K.  An  equation  similar  to  (2)  but 
applicable  here  is 

w  =  5.0  x  10"5  +  28  q  (13) 

u  u 

The  900  mb  surface  was  assumed  to  be  a  black  body.  Eq.  (3)  was  then 
evaluated  at  300,  500,  700  and  900  mb,  and  cooling  rates  were  com¬ 
puted  at  400,  600  and  800  mb.  Existence  of  cloud  at  these  levels 
was  deduced  using  values  of  r(p)  (see  (12))  obtained  in  the  same  manner 
as  those  in  Table  1 .  Surface  observations  of  clouds  were  also  utilized. 
In  other  respects,  the  procedures  of  Section  3  were  followed. 

Results  are  presented  in  Figs.  6-9.  The  reader  is  referred  to  the 
paper  cited  above  (Danard  1964)  for  synoptic  charts  and  other  details 
of  the  case  studied. 

Fig.  6  shows  the  net  flux  (positive  upward)  at  300  mb.  The  top 

of  the  principal  cloud  deck  is  clearly  delineated  by  relatively  low 

values  of  flux.  Its  lateral  boundary  is  given  approximately  by  the 

5  -1  -2 

isopleth  F  =  1.75  x  10  ergs  sec  cm 

Figs.  7-9  show  the  rate  of  loss  of  energy  due  to  long-wave 
radiation  from  water  vapor.  By  far  the  greatest  loss  occurs  in  the  upper 
portion  of  the  cloud  mass.  The  large  values  in  the  Oklahoma-Kansas 
region  at  800  mb  (Fig.  9)  indicate  strong  cooling  at  the  top  of  strato- 
cumulus  cloud  with  clear  skies  above.  Zero  values  mean  that  the 


11 


point  in  question  is  in  cloud.  Negative  values  indicate  warming  be¬ 
neath  cloud  or  in  the  lower  portion  of  the  cloud  mass .  For  comparison 

-1  2  -1  -1 

with  the  results  of  Section  4,  1  C  day  =  1.17  x  10  ergs  gm  sec 
Figs.  7-9  are  to  be  compared  with  Fig.  8  of  Danard  (1964).  The 
latter  depicts  the  rate  of  release  of  latent,  heat  at  400,  600  and  800  mb. 

This  was  largest  in  the  Indiana- Pennsylvania  area  and  amounted  to 

3-1-1  3 

about  3  x  10  ergs  gm  sec  at  400  and  800  mb  and  about  5x10 

ergs  gm  ^sec  1  at  600  mb.  Thus  in  the  case  studied,  the  long-wave 
cooling  at  cloud  top  was  smaller  by  a  factor  of  about  five  than  the 
maximum  rate  of  release  of  latent  heat  in  the  upper  troposphere.  It 
ma>  be  remarked  that  this  case  was  one  of  heavy  precipitation  and  in 
less  intense  cyclones,  the  ratio  would  be  smaller.  However,  In  the 
middle  and  lower  troposphere,  radiative  loss  was  generally  negligible 
in  the  region  of  precipitation  compared  to  release  of  latent  heat. 

The  long-wave  cooling  at  cloud  top  amounts  to  5.4  C  day  1  in 
the  upper  200  mb  of  the  cloud  mass  (Fig.  7),  This  latter  figure  agrees 
qualitatively  with  the  estimate  of  Moller  (1951).  The  cooling  would 
tend  to  be  offset  by  absorption  of  solar  radiation  at  cloud  top.  How¬ 
ever,  many  workers  (e.g.,  Moller  (1951)  or  London  (1957))  regard  this 
effect  as  small.  The  long-wave  cooling  at  cloud  top  greatly  exceeds 
that  in  the  cloud-free  air.  Thus  an  accurate  computation  of  the  latter 
may  be  relatively  unessential  in  the  study  of  cyclogenesis  over  short 
time  intervals  (say  48  hrs  or  less). 

The  temperature  difference  between  800  and  400  mb  at  the  location 
of  the  sea-level  low  (southern  Illinois)  is  31  C.  By  comparison,  if 
the  lapse-rate  were  pseudo-adiabatic,  the  value  would  be  34  C. 
Neglecting  effects  of  vertical  motion  and  differential  advection,  and 
assuming  that  the  pattern  of  computed  radiative  cooling  moves  with 
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the  cyclone,  the  800-400  mb  temperature  difference  at  the  surface 
low  would  be  36  C  after  24  hrs.  This  represents  a  conditionally  un¬ 
stable  state  which  would  presumably  be  alleviated  by  convective 
overturning. 


6.  Concluding  remarks 

A  possible  modification  to  the  treatment  of  clouds  would  be  to 
allow  for  fractional  amounts  in  each  layer.  However,  this  would 
complicate  the  calculations  unless  the  number  of  layers  in  the  vertical 
were  kept  small. 

Short-wave  absorption  by  water  vapor  is  appreciable  in  the  lower 
troposphere  at  low  latitudes  (see  Section  4).  The  average  heating  due 
to  this  source  could  be  included  if  desired. 

The  question  of  the  influence  of  long-wave  radiation  on  the 
evolution  of  cyclones  is  probably  best  answered  by  comparative  num¬ 
erical  integrations.  Cooling  and  destabilization  of  the  upper  part  of 
the  cloud  mass  are  to  be  expected.  Destabilization  and  subsequent 
increase  in  precipitation  may  be  important  in  the  dynamics  of  de¬ 
velopment. 
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Legends  for  Figures 

Fig.  1.  Schematic  diagram  showing  levels  referred  to  in  Eq.  (3). 

Fig.  2.  Schematic  diagram  showing  levels  at  which  terms  in  Eqs.  (3) 
and  (9)  are  computed  in  Section  3. 

Fig.  3.  Configurations  of  cloud  in  the  model  atmosphere  of  Section  3. 
Hatching  denotes  cloud. 

Fig.  4.  Annual  average  cooling  rates  f'ra  cloudless  atmosphere 
associated  with  long-wave  radiation  from  water  vapor. 

Units:  deg  C  day  1 . 

Fig.  5.  Difference  between  cooling  rates  of  Fig.  4  and  those  of 
Manabe  and  Moller  (1961),  who  include  both  long-  and 
short-wave  radiation.  Units:  deg  C  day  *. 

Fig.  6.  Net  upward  flux  at  300  mb.  1200  GMT  21  January  1959, 

tt  ,n5  -1  -2 

Units:  10  ergs  sec  cm 

Fig.  7.  Net  loss  of  radiative  energy  per  unit  mass  at  400  mb. 

2  - 1  - 1 

1200  GMT  21  January  1959 .  Units:  10  ergs  gm  sec 
Fig.  8.  As  in  Fig.  7  but  for  600  mb. 

Fig.  9.  As  in  Fig.  7  but  for  800  mb. 
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